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We study Cerenkov radiation from moving straight strings which glisse with respect to each other 
in such a way that the projected intersection point moves faster than light. To calculate this effect 
we develop classical perturbation theory for the system of Nambu-Goto strings interacting with 
dilaton, two-form and gravity. In the first order one encounters divergent self-action terms which are 
eliminated by classical renormalization of the string tension. Cerenkov radiation arises in the second 
order. It is generated by an effective source which contains contributions localized on the strings 
world-sheets and bulk contributions quadratic in the first order fields. In the ultra-relativistic limit 
radiation exhibits angular peaking on the Cerenkov cone in the forward direction of the fast string in 
the rest frame of another. The radiation spectrum then extends up to high frequencies proportional 
to square of the Lorentz-factor of the relative velocity. Gravitational radiation is absent since the 
1+2 space-time transverse to the straight string does not allow for gravitons. A rough estimate of 
the Cerenkov radiation in the cosmological cosmic strings network is presented. 

PACS numbers: 11.27.+d, 98.80. Cq, 98.80.-k, 95.30. Sf 



I. INTRODUCTION 

During past few years the hypothesis of cosmic strings received new impetus from superstring theory. Although 
perturbative superstrings in ten dimensions are too heavy to be admitted as cosmic strings, it was reaHzed that there 
are various possibiHties for geometry of string compactifications to accomodate four-dimensional strings with much 
lower tensions. Copious production of cosmic strings is typical [H, 01 for the brane inflation scenario [sl, in which the 
period of inflation is associated with the coUision of branes. This scenario provides an acceptable model of inflation 
and predicts creation of cosmic superstrings consistent with the current CMB data. These strings typically have 
lower tensions than the usual GUT cosmic strings and thus they are not the main players in the formation of cosmic 
structures, but their observational signatures could provide a direct conflrmation of the string theory. This stimulated 
a detailed study of creation and evolution of the cosmic superstring network within the KKLMMT model [3], racetrack 
inflation Q and other particular scenarios 0, 0, (for recent review see [£, Jfll, [HI, One particular feature 

of the warped IIB compactiflcations involved in these considerations is prediction of two types of cosmic strings: 
F-strings (fundamental) and D-strings (Dirichlet) with different tensions /if, (3, [H, [il, F and D strings 
may also form the so-called p, q— composites [18] and provide for triple junctions with new exotic observational 
predictions. Typical values of the dimensionless parameter G/U are in the interval (10~^^ — 10"^). Another new feature 
is that cosmic superstrings generically have lower reconnection probabilities lying in the range P ~ (10"'^ — 1). This 
changes cosmological evolution of the string network j^S, lU, [23, [H, [131 leading in particular to enhancement of the 
fraction of straight strings. 

Both the string network evolution and their possible observational signatures crucially depend on radiation processes. 
In fact, it was re cog nized long ago that oscillating loops of cosmic strings generate large output of gravitational 
waves fiF, '2^, "2?, '28] at the level accessible for current and future detectors. Global strings produce massless axions 
[28 . 29 . 30 , 3L 32., ^] which become massive at a later stage of expansion, creating an observational constraint on 
the axion mass. In the models containing the massless dilaton (like cosmic superstrings), the dilaton radiation from 
strings may also constrain the string tension parameter Gfi [s^. Issl. [s^. 

The main mechanisms of radiation in the string network which have been explored so far were radiation from 
smoothly oscillating string loops and radiation from kinks and cusps formed on them. It was tacitly assumed that 
straight unexcited long strings do not radiate. Meanwhile, interaction of long strings via massless flelds gives rise to 
another radiation mechanism of Cerenkov nature. When two straight non-excited Nambu-Goto strings interacting at 
a distance via dilaton, two-form and gravity glisse with respect to each other, they get deformed in the vicinity of the 
point of their minimal separation. This point can propagate with a faster-than Hght velocity, provided the incHnation 
angle between the strings is sufficiently small. (For the strictly parallel strings this velocity is infinite for any relative 
string speed.) In this case, the propagating deformation together with the associated field tensions become the source 
of Cerenkov radiation, which is the effect of the second order in interaction of the string with massless fields. This 
mechanism was suggested in [s^l for gravitationally interacting strings, but it turned out that, although the effect is 
kinematically allowed, the corresponding amplitude is zero on the mass-shell of the graviton. The reason is that the 
1-1-2 space-time orthogonal to the straight string does note allow for gravitons. (However, gravitons are produced at 
quantum level in the strings recombination and annihilation [s^l-) But Cerenkov mechanism works for other possible 
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string massless excitations leading, e.g., to electromagnetic radiation from superconducting strings [39| or emission of 
axions. The latter efect was recently studied in detail in the fiat space-time |40|- It was found that, although being of 
the second order in the axion coupling constant, it still gives a considerable contribution into the total cosmological 
axion production by the string network. 

Cerenkov radiation from straight strings is similar to bremsstrahlung of point charges in electrodynamics. Moreover, 
the system of parallel strings interacting via a two-form field in D space-time dimensions is exactly equivalent to 
the system of point charges interacting with the vector field in (D-l)-dimensional space-time. A distinction of the 
Cerenkov mechanism of dissipation in the cosmic string network from conventional radiation via formation of loops 
can be understood as follows. Formation of loops from initially disconnected straight strings is effected via direct 
contact interaction of intersecting strings. In our case the long-range interaction of strings via massless fields which 
can be potentially radiated underlies the formation of the superluminal radiation source. Cerenkov radiation is a 
higher order effect and thus potentially smaller. But it works for a wider set of initial data in the string network 
(non-zero impact parameters of colliding strings), than that corresponding to intercommuting strings. Also, as we will 
show here, Cerenkov radiation is strongly enhanced in the case of relativistic velocities. So it can still be non-small in 
the cosmological setting. The detailed cosmological applications remain outside the scope of the present paper, but 
we give the rough estimates of the dilaton and two-form field backgrounds generated via Cerenkov mechanism in the 
evolving string network. 

In this paper we consider Cerenkov radiation of moving straight strings interacting with three massless fields: dilaton, 
antisymmetric second rank tensor (NS-NS or RR two-forms) and gravity. To facilitate construction of the solution 
of the coupled system of the Nambu-Goto and the field equations in the second order we introduce a diagrammatic 
representation similar to Feynmann graphs. In the first order approximation one encounters divergencies which are 
cured by classical renormaHzation of the string tension. Then the effective sources of the dilaton and two-form radiation 
are constructed and the radiation rates are calculated. 

The plan of the paper is as follows. In Sec. 2 we give the general setting of the problem and introduce graphic 
representation for classical vertices corresponding to interactions of strings with dilaton, two-form and gravity and to 
interactions between the fields. Sec. 3 contains formulation of the perturbation approach and an explicit construction of 
the second order equations. In Sec. 4 we reproduce in our framework the main features of classical renormaHzation for 
strings (in the lowest order of perturbation theory) : elimination of scalar and two-form divergences via renormaHzation 
of the string tension, and the absence of gravitational divergency. In Sec. 5 the first order deformations of the string 
world-sheets due to interaction via the dilaton, two-form and Hnearized gravity are considered. The main calculation 
is presented in Sec. 6 where we construct radiation generating currents by computing the source terms in the wave 
equations of the second order. These include the world-sheet local terms and the bulk terms which correspond to 
multiple second order graphs including all relevant elementary vertices. Radiation rates for the dilaton and the two- 
form fields are computed in the Sec. 7 and analyzed in detail in the case of ultrarelativistic velocities. FinaHy, in Sec. 8 
we present rough cosmological estimates and summarize our results in Sec. 9. 



II. ACTION AND EQUATIONS OF MOTION 

We consider a pair n = 1, 2 of straight Nambu-Goto strings described by the world-sheets 

x^=X^{a^), M = 0,l,2,3, a, = (T,a), a = 0, 1. (1) 

Strings interact via the gravitational 5^1/ = + hf_i^, the dilaton (t>{x) and the two-form (axion) B^^{x). Using the 
Polyakov form for the string action, we present the total action in the form: 

+ j i^2d^cj,d,^g^,+ ^-H^,pH^'^Pe-^'^0'-^^^^gd^x, (2) 

where the sum is taken over n — 1,2 (the corresponding index is omitted for brevity). Our signature choice is mostly 
minus for the space-time metric, and (-1-,—) for the world-sheets. The Levi-Civita symbol is e*^^ — — e^" — 1, jab is 
the metric on the world-sheet, and the lower Latin indices mean partial derivatives with respect to the world-sheet 
coordinates Xj^ = daX^^ . The action contains four parameters: the string tension ^, the Newton constant G, the 
dilaton coupling a and the two-form coupling /. The antisymmetric three-form field strength is defined as 



(3) 
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Variation of the action ([2]) with respect to X^^ leads to the equations of motion for strings 



The coordinate derivatives on the right hand side apply to the metric, dilaton and two-form fields, their values are 
taken on the world sheet. The derivatives with respect to the world-sheet coordinates cr° on the left hand side apply 
both to the world-sheet quantities X'''{a"'), jab and to the metric, dilaton and two-form fields, e.g., for the metric 

c'a.g^i. = X^dxg^u, (5) 

and similarly for 0, B^^. In this formalism "fab is an independent variable; variation of the action with respect to "fab 
gives the constraint equation 

{X^:X^, - ]^"fabl"'Xi^X'^^)g^, = 0, (6) 

whose solution defines "fab as the induced metric on the world-sheet: 

"lab = X^X'^g^,. (7) 
Consider now the field equations. Variation over gives the dilaton equation 

(s^'-a.^V^) + ^ij2e-4a0 ^-T.tJ ^l^^bg^.^r'V^e^'^^SHx - X{a, r))d^a, (8) 

where the sum in the source term is taken over the contribution of two stings. The equation for the two-form field 
reads: 

{H'^'^^e-'^'f'V^) =-J2^7:fJ X:X>;e'^H\x - X{a, T))d'a. (9) 
We also have the Bianchi identity 

V[^i?a/3T,] = 0, (10) 

where alternation over indices has to be performed and the derivative can be equivalently treated as a covariant 
derivation with respect to g^i, or a partial derivative. 
Finally, for the metric we have the Einstein equations 

R^.u - ^g,.,^R = 87rG(r,„.+ 1^^+ T^^), (11) 



where the source terms read 



St 



Y: M / X,^X,-7°^7^e^"^ ^ dV, (12) 

T'^'' = 4 fa^^a^ - Igf^-^i^^f) , (13) 



2- 

The total system of equations consists of two equations for stings of the type ^ and the field equations |[8l O [TT|) . It 
describes classically in a self-consistent way motion of the strings with account for their interaction via the dilaton, 
two-form and gravitational fields as well as evolution of the generated fields. This system can be solved iteratively in 
terms of the coupling constants a, f, G. To facilitate this construction it is convenient to introduce a classical analog 
of the Feynmann graphs, denoting string by bald solid line, the dilaton by this solid line, the two-form by dashed and 
graviton by wavy lines (Fig. 1). Simple analysis of the equations reveals that we have vertices involving linear and 
non-linear string-field interactions, as well as three-leg and multi-leg interactions of fields between themselves. 
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Phc. 1: Vertices associated with the action graphs a,b,c are string-dilaton, string-two-form and string-graviton vertices of 
the lowest order; graph d depicts higher-order string-dilaton-graviton vertices (for all k > 1, k' > 1); graphs e,f,g are lowest 
order field interactions present in the action ([2|), and graphs e,f,g are multi-graviton vertices accompanying the lowest order 
ones. 

III. PERTURBATION THEORY 

We construct an iterative solution of the coupled string-field system expanding the string world-sheet mapping 
functions X,j^{(j°-), n — 1,2, and the field variables in terms of the coupling constants a, f, G: 

- i + l + 

1 2 
B — B^V^ B^u ^ . ■ ■ , 

1 2 

hf,y = ht_lV+ hpLU + (15) 

Here the expansions of start from the zero order, while those of the field variables — from the first order terms, 
that is we assume that there are no background dilaton, two-form or gravitational fields. Zero order mapping functions 
describe the straight infinite uniformly moving strings 

h = d^n + <r (16) 

Here is the unit spacelike constant vectors along the strings, and are the unit timelike constant vectors — four- 
velocities of the strings. The corresponding three-dimensional velocities are orthogonal to the strings. The constant 
vectors can be regarded as impact parameters for two strings with respect to the chosen frame. In the zero order 
the space-time metric is fiat, and the corresponding world-sheet metrics are also Minkowskian ryafc = diag(l,— 1) in 
view of the normalization assumed 

(!]„£„) = Ty^t^SjJSj'j ^ -1, {UnUn) = 1, (SnUn) = 0. (17) 

It is convenient the choose the Lorentz frame in which the first string is at rest and is stretched along the z-axis: 

< = [1,0,0,0], = [0,0,0,1]. (18) 
The second string is assumed to move in the plane , with the velocity v orthogonal to the string itself: 

^2 = 7 [1, 0, — w cosQ!, t; sin a], E2 = [0' siK^: cos a], (19) 

where 7 = (1 — u^)~^/^. In such a configuration the strings never intersect each other remaining always in the parallel 
planes. Apart from the orthogonality conditions (flTl) . four other scalar products are 



(wi,M2)=7, (El, S2) = - cosa, (mi,S2)=0, (u2, Si) = -U7sina, 



(20) 
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note that ui and S2 are orthogonal. We also choose both impact parameters to be orthogonal to and and 
aligned with the axis , the distance between the planes being d = d2 — di. The angle of inclination a of the second 
string with respect to the first one can be written in a Lorentz-invariant form 

a = arccos(-i;iI]2). (21) 

Similarly, the invariant expression for the relative velocity of the strings is 

V = (1- (uiU2)-')i (22) 

With this parametrization of the unperturbed world-sheets, the projected intersection point (the point of the minimal 
separation between the strings) moves with the velocity 



sma 



1 

I ._i / (uiM2)^--i y 



along the x^-axis. This motion is not associated with propagation of any signal, so the velocity Vp may be arbitrary, 
in particular, superluminal Vp > 1. The case of parallel strings corresponds to Vp = 00. 

Note that the above introduced parameters l|22|2ip are not invariant under reparameterizations of the world-sheets 
[4o| . The quantity which is invariant under the volume preserving reparameterizations is 

K = det (X'^^X^^ijf,^) = 7 cos a. (24) 

The superluminal regime corresponds to k > 1. 

The expansions (|15p are substituted into the system of equations |[4l [Sj [9l llip which has to be solved iteratively. 
The zero order differential operator in the dilaton equation (jlj is the fiat-space D'Alembert operator □ = —if'^d^dy. 
Similarly, choosing the Lorentz gauge for the two-from and the metric perturbations 



d^B^^^Q, d^V"^^, ^M.^ = /iM-_^M-/i, (25) 



where h = h^, we get the linear D'Alembert equations for the first order two-form and the gravitational field as well. 
Due to Hnearity of the field equations, the first order dilaton, two-form and metric perturbations can be presented as 
the sums of the separate contributions due to two strings 



(26) 



Here each term with n — 1,2 satisfies the individual D'Alembert equation with the source labelled by the same index 
n: 



1 


1 1 




= 01+ <fe, 


k" 


= sr+^r 







□ i„ = 47r J„, (27) 



□ = 47r jr, (28) 

□ = 47r . (29) 

The coupling constants are included into the source terms, while zero indices in the sources indicate that they are 
computed using the zero order approximations for the strings mapping functions. The source terms thus read 



Jr. = ^ I (x~ Xn (t,(t)^ dV, (30) 



h ^ U K'S'ix-Xi T^a)] d'a, (31) 



2 

4GAi 



„ J Wii^S^ l^x~ Xn (r, a)^ dV, (32) 
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where the following antisymmetric and symmetric tensors are introduced 



,0 
ab T^/j ^1/ 



bn 



t^n"^ ~ Xan ^bn = ""n^n ~ ^ri^nJ Ul^'^Tjuu ~ Un — 2, 

Now consider the perturbations of the strings world-sheets induced metrics (index n is omitted) 



(33) 
(34) 
(35) 



(36) 



Jab = 11ab+ lab, 

where in the general gauge the first order correction reads 

1 00^1 "^1/ 

lab -^b hp.v 4^ 2 JCbVl^Li^- 

However, to be able to disentangle the higher order perturbed equations one has to get rid of the second term in l(37|) 
1 

proportional to X" ■ Using the space-time and the world-sheet diffeomorphism invariance, we can impose the gauge 
condition 



(37) 



1 ^ 
^a XbVfJ.i' = 0: 



(38) 



in which case the perturbed induced metric will contain only the zero order string mapping functions: 

lab =h h 1= U^-" (39) 



The first order perturbations of the mapping functions describing deformations of the flat world-sheets satisfy 
the following equations following from the Eq.Q: 



hf^u + \yhu 7 ) rf''-T]f,^ 1 



yah 



where raising of the world-sheet indices is performed by the Minkowski metric: 

=icdv""v'"', i^icdv"^- 

Differentiating the field variables along the world-sheet according to the rules 

10,1 1 0,1 1 0,1 

da (f>^X^dx (t), da Bt,u --Xadx Bf,^, h^^ =Xadx ht,u, 



we can rewrite the above equation as 



^^rj^^dadbX'-^Fl^^^+F^^^+Fl^,.^, 



where the forces due to dilaton, two-form and graviton are introduced: 



1 V 



F^B) = '^^fy' Kx, 



hXr... i U^^W"^ - -U^^'U^^ ]~U,xi /l^"'^ 



(40) 
(41) 
(42) 
(43) 

(44) 

(45) 
(46) 



In these equations the indices labelling the strings are not shown, but it is understood that for each string we have to 
take into account at the right hand side both the self-force terms arising form the fields due to the same string, and 
the mutual interaction terms coming from the partner string. 
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Consider now the second order field equations. Expanding the Eqs. l(8l)- lfTT|) to the next order and imposing again 
the Lorentz gauge for the two-form and the Unearized gravity one obtains the D'Alembert equations for the second 
order fields with the source terms involving contributions due to deformations of the strings world-sheets (local terms) 
as well as the quadratic combinations of the first order fields (bulk terms) . The dilaton equation with account for ([39|l 
will read 



- 4 X'^dn 



x~ X (r, cr) 



— H ' Hfj,i/X — dfi 



^V^"" h ] du 



(47) 



where the partial derivative operator acts on the delta function. Note that the first order terms in the first line are 
multiplied by the dilaton coupling constant and thus give the same order quantities as the products of two first order 
field quantities. 

The equation for the second order two-form field is: 



□ 1^'^ ^ 47r = 



X- X (r, cr) 



dx 



1 1 1 
- h -4a (/} 



H 



XflV 



(48) 



where alternation over indices includes the factor 1/2. 

The situation is slightly more complicated for the graviton. To obtain an equation for the second order gravitational 
perturbation one has to take into account the quadratic terms in the Einstein equations. In the gauge 

one has the following expansion of the Einstein tensor up to the second order: 



1 1 11 

where -0^,^ contains the first and the second order quantities, and the quadratic term reads: 



l-^l^.A'2h"^^^^xh^0 - dxh^pdph''^ 



-dxh^pd^h'^'') 



(49) 



(50) 



Extracting the second order terms we obtain from the Eq. (fTT|) : 



where the source term reads: 



= f E/ (2l(^i.)a-i'5A)5n^-l(r,cT))rfV + 

+ f E J [{U^,,U^^ ~2U^U:) hxr +U^,,{4a'^ - h)]s''{x^ X iT,a))d^a + 

1 11 1^111 1 ^ 1 1 1 > 

+ G ^S^,4hxT) + 2df, (t)d^(j) -ri^^^dx <l> cj) +- H^xr H^^^ - —v^u H^xr H"^'' . 



(51) 



(52) 



All source terms in the second order field equations have similar structure. Note the presence of the derivatives 
from delta functions in the strings world-sheet contributions. 

It is worth noting that the second order dilaton, two-form and gravity fields already can not be presented as a 
sum of terms due to separate strings, but rather look as generated by the collective sources. These sources contain 
contributions not only from the perturbed world-sheets, but also the bulk contributions which are not associated 
with separate strings. As we will see later, these sources may have superluminal nature, in which case the Cerenkov 
radiation will appear. 
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a b c 



Phc. 2: Graphs describing self-interaction due to dilaton (a), two-form (6) and linearized gravity (c). The contributions from a 
and h are divergent and have different signs. The contribution of c is zero. 



IV. RENORMALIZATION 



The action of the proper fields upon the source string is described by the self-action terms in the Eq. l(43|) 
corresponding to the graphs shown on Fig. 2. the dilaton and the two-form lead to divergences while the contribution 
from the lineariz ed g ravity vanishes. Divergent terms can be absorbed by renormalization of the string tension 
parameter ^ (isl . HJ We consider renormalization in the first order of the perturbation theory. Linearizing the 
string part of the action ([2]), one can split it into the sum: 

S — Sst + + Sb + Sh, (53) 

where Sst is Polyakov action with the bare tension parameter: 

Sst = -y / xiixi^v^..v^'d^<J, (54) 

and three other terms describe interaction with the dilaton, the two-form and the linearized gravity: 

= -Moa J X^X^i^^.r^'^'' So, (55) 

Sb = -2TTf J B^, X^X'^^e-" d^a, (56) 

Sn ^ j X^^X^.^^^h^, dV. (57) 

Since we are working in the lowest order of the perturbation theory, the mapping functions here are quantities of 
zero order. In obtaining the last formula we used the following expansion: 

7"''%/^ = + V Q^c^^'' -X^'''X''^^+■■■ (58) 
Consider the first order dilaton field on the world-sheet of the source string 

a^o /■ e*«-('*"-^"(^'^_»%u)%E2 4 

2 "'Q.J V^"j 



' ' 8tt^ J q 

where qu — = g^S'' are the fiat-space scalar products, and — q^^q^ ■ Due to delta-functions, the scalar 

product in the exponential is constant: 

g^X''(r, cr) = {qd) + {qu)T + (gE)cr = [qd) = const, (60) 

so the integrand does not depend on a and r. The integral diverges as 

5{qu)6{q^) ^^f^-dq^ (g^^ 

9 Jo qi- 

where we used the frame (fTSl and introduced polar coordinates in the 2-plane orthogonal to and S'' : 

q^ = (qu)^ - (q^)^ - ql ~ ql qj + qj ^ ql. (62) 

The integral logarithmically diverges at both ends. With the infrared (IR) an ultraviolet (UV) cut-off parameters 
g™'", qf'^^ one can write: 

„max 

/==27rln^^. (63) 
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Substituting this into the Eq. (|59|) and further into (|55|) we find the regularized dilaton part of the action: 

,2^,2 

In 



ii^n'^ „max /• 

= ^^-'p^J X^^XM'd'-- (64) 



Since the functional is the same as the bare action functional Sgt^ one can absorb divergencies by renormalization of 
Mo- 

Similarly, the first order two-form field on the world-sheet of the source string reads: 

(^''"^ = ^7 -q2 1- (6^) 

This integral also diverges as l(63|) . In view of the relation 

yl^^^yX = (66) 

one can see that the action (|56|) also have the functional form of l(54|) . namely 

5^^^ = -2nf In / XiiXi;r^,^ri^''d^a. (67) 

Finally for the graviton (r, cr) on the world-sheet we have the divergent integral 



where Wf^^ id given by the Eq. l(35|) . However, substituting this into l(57|) one obtains zero in view of the identity 

Wf^.U'^" = 0. (69) 

Therefore, gravitational interaction of the strings does not lead to classical divergences in the lowest order of the 
perturbation theory. This result is conformal with previous results [ill, obtained with different tools. 

Collecting the above formulas, we see that to remove self- interaction divergences one has merely to replace the 
tension parameter in the action ((54| as follows: 

Divergences due to the dilaton and the two-form have opposite signs. This refiects the fact that the scalar interaction 
is attractive, while interaction via the two-form is repulsive. If the Bogomolny-Prasad-Sommerfield (BPS) relation 
between the dilaton and the two-form caoupHngs is satisfied 

4^=4./^ (71) 

the divergent terms mutually cancel, and there is no renormalization at all, fi — no (for earlier work on this subject see 
refs. [41, 42, 43, 44, 45, 46]. Note that our dilaton coupHng constant has dimension of length, the usual dimensionless 
constant a (quantity of the order of unity) is related to it as 

= Ga^. (72) 

It has to be noted that renormaHzation can be performed in a simple way only at the linearized level. When all 
non-linearities are taken into account, classical renormaHzability of the bosonic string theory interacting with gravity 
is lost. In this paper we will be restricted by the second order calculation of radiation which involves only on-shell 
second order quantities. These are unaffected by the higher order renormalization effects. So, in what follows, we will 
use the value for the string tension renormaHzed in the first-order, and thus omit all self-interaction terms. 
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Phc. 3: Deformations of the strings world-sheets due to interactions via the dilaton (a), two-form (b) and linearized gravity (c) 
n = 1,2, n' = 1,2, n / n'. 



PERTURBATIONS OF THE STRING WORLD-SHEETS 



Now consider the first order perturbations of the mapping functions X^(t, cr), n = 1,2 caused by mutual 
interactions. We have to substitute into each string equation of motion ^ the first order fields generated by another 
string. It is convenient to use the following index convention: n = 1,2, n' = 1,2, n' . The total perturbation thus 
splits into three terms due to dilaton, two-from and graviton exchange: 



I ip B h 



^j:, (73) 

as shown on Fig. 3. Let us start with dilaton exchange contribution (Fig. 3a). The corresponding world-sheet 
perturbation X^(t, a) is the solution of the two-dimensional D'Alembert equation (following from the Eqs. I|43|44p ): 



-2a 



[u^^d, - a^) (x) 



J a;=X„(T,(T) 

where the dilaton field generated by the string n' is taken on the world-sheet of the string n: 



da 



x=X-a{T,(y) S,TT''l 



^S{qUn')S{qT,n' 



g2 -I- 2ieq° 



qf, d'^q. 



(74) 



(75) 



Substituting l|75l) into l(74|) one can obtain the solution dividing the right hand side by the two-dimensional D'Alembert 
operator as follows: 



/i [ A„, [g^ + S^J(gS„) - <(gM„)]e-'9^^(r-) 



(27r) 



q^iiqun 



where 



An' =e''<'6iqUn')SiqJ:n'). 



(76) 



(77) 



Note, that delta-functions in the integrand have support outside both the Hght cone = and the surface {qu,, 
{q^n)"^ (except for the trivial point q^ = 0), so the integral is finite. 

Consider now the two-form interaction (Fig. 3b). We have to solve the equation 



[di - di) h = 



2:=X„(T,cr) 



where the three-form is 



x=Xn{T,a) 



27ri 



g2 + 2ieqO 



(78) 



(79) 



where curly brackets {} denote the cyclic permutation of indices. Again, dividing by the operator (5^ — 9^) one 
obtains: 



Similarly, the gravitational contribution is described by the equation: 



(52 - dl) x^ = ^u. 



9^ L/3 -29o kp - Ur [d, L/3 -2dp Lu 



x=X,i{t,<t) 



(80) 



(81) 
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where the variation hf_ii, is generated by the partner string n': 



1 



g2 + 2ieq° 



d^q. 



(82) 



Solving this equation one finds: 



-d\. 



(83) 



It can be checked that the gauge condition l|38p imposed in the beginning of the calculation holds indeed for each 
of the three separate contributions to the perturbed mapping functions. 

VI. EFFECTIVE SOURCES OF RADIATION 

1 1 1 

The first order fields (j), B^" , /i^" do not contain the radiative parts. Consider, e.g., the Fourrier-transform of the 
dilaton: 



4>{k) = / (j){x)e"'''d^x 



(84) 



The retarded and advanced solutions of the first order wave equation l(27|) will read: 



(/)„ (fc) — 2tt^ 



± 2iek° 



(85) 



The radiative part 



k^ik) = K kik)- kik)) = -in^^la6iku^)S{kJ:„)S{e)e''''i" 



(86) 



is the distribution having support only at the trivial point = in the momentum space. Thus, to investigate 
radiation, we have to pass to the second order of the perturbation theory. The problem reduces to the construction 
of the source terms in the wave equations of the second order. 



Dilaton 



Consider the second order dilaton equation l|47p in more details. The current at the right hand side contains the 
contributions localized on the string world-sheets (the upper line in l(47l) ) and the bulk contributions coming from 
the products of the first order dilaton, two-form and graviton fields (the lower line). The former contains the sum 
over the strings which can be understood as follows. One has to take the perturbations of the world-sheet mapping 
functions X^^ for each string n = 1,2 due to first order field generated by the partner string n = 2,1 respectively. 
These contributions are depicted by the graphs a„, 6„, c„, n— 1, 2 on the Fig. 4. The external dilaton leg corresponds 
to the emitted dilaton with the momentum k^ (in our classical treatment to the Fourrier-transform of the current) , 
thus the terms in the sum in l|47p ) with n = 1 are given by the graphs ai,bi,Ci and those with n = 2 — by the 
graphs a2, 62, C2. Other upper line terms must be treated in a similar way: one has to take i/^'' for each n = 1, 2 and 
multiply by the graviton and dilaton perturbations caused by the other string n = 2,1. These are depicted by the 
graphs dn, e„, n = 1,2. In this way the self-action of the fields upon the string will be excluded. 

On the contrary, the terms in the lower line of the Eq. (|47|) are non-local and not pairwise. Here the self-action 
terms also have to be excluded, so we take in the quadratic combinations only the products of the first order fields 
generated by different strings. 

Consider first the contributions described by the graphs oi, bi, ci (dilaton emission from the first string line). The 
corresponding current reads: 



St 

Jl 



(x) 



"8n 



J X'^d^S^[x~ Xi{t,<j) 



d^a. 



(87) 
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2 




f g h 

Phc. 4: The diagrams contributing to dilaton radiation in the second order of the perturbation theory: ai, 6i, ci correspond to 
deformation of the first string and 02, 62, C2 - to deformation of the second one; graphs d, e stand for contact terms. Diagram / 
corresponds to the product of (first order) two-form fields generated by the first and the second string, graphs g, h — to mixed 
graviton-dilaton contributions. 



where the perturbation of the mapping function of the first string Xi is caused by the first order dilaton, two-form 
and graviton fields generated by the second string. Substituting the corresponding terms into (|87|) we obtain after 
some rearrangements: 



The vector is the sum of three terms according to the decomposition described above: 

- + '-j^rr + ^zr, (89) 

where the dilaton exchange contribution is 



7^5" = g" + E5'(<zSi)-<(gui), (90) 

the two-form contribution is 

r/' = [(?/iU2)(SiE2) - (SiW2)(«lS2)] + [(gui)(M2Si) - {^iq){uiU2)] + 

+ u^[(uiS2)(Si9)-(gTti)(EiI]2)], (91) 

and the graviton contribution is 

= [{uiU2f + (S1S2)' - ("1^2)' - {U2^lf - 2] + [uliqm) - E5'(gEi)] • 

[{uiU2f - (EiS2)2 - (wiE2)2 + (7/2Ei)2 + 2] - 2 K(gEi) + Si'(gui)] [{uiU2)iu2^i) - (SiS2)(7iiS2)] - 
- 2ut^ [{uiU2){qui) - (EiU2)(gEi)] + 2E^ [(wiE2)(gui) - (EiE2)(gSi)] . (92) 

Consider the Fourier-transform 

J{k) = / J{x)e'''^d'^x (93) 



with fc^ = (on-shell condition for the massless particle) . Substituting the above expressions and integrating over the 
world-sheet of the first string one obtains two more delta-functions: 

g,(fc_,)(„i^_Si<T)^^^^ = {2nf S[{k - q)ui]5[{k - g)Ei], (94) 
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so totally we will have the product of four delta-functions in the integrand: 

Aiiq, k) - S{qui)5{q^i)d[ik - q)u2]S[ik - g)E2]. (95) 

Now consider the contribution J2{k) coming from the second string (the graphs a2,62,C2 on Fig. 4). Obviously it can 
be obtained from the previous result by interchanging indices 1^2 elsewhere. In this case we will get the product of 
the S -functions in the form 



A2(g, k) = S{qu2)S{qJ:2)S[{k - q)ui]S[ik - g)Si 



(96) 



It is convenient to cast this second integral into the same form as the previous one. For this it is sufficient to shift the 
integration variable in J2(k) as follows: q^ k^^ — q^^ . Since ^.l{q^ k) = A2{k — q, k), we will get the same product of 
the 5-functions, so finally we can present the total contribution of the first six graphs as follows 



j{k)=j, {k)+j2 {k)= I n(g,fc) 



e i (g) ^ 62 (fc - q ) 
g2 {k ~ qy 



(97) 



where 



n(g, k) = e'('iifc+'?('i=-'ii))Ai((j, fc), 



(98) 



and 



St 

81 (<?) 



ikDi)+T:f{kYi) + Gii^ikZi) 



(99) 



To get the function 82 from Qi one merely has to interchange the indices 1, 2 labelling vectors in the scalar products 
changing simultaneously g'' ^ (fc — g)^. 

The contributions of the next four graphs di,d2, 61,62 are computed similarly. The resulting "contact" term in the 
source can be presented again in the form Ip7|) with equal contributions from two strings 



61 (g) =62 (g) = ^ + [("1^2)' + (SiS2)2 - («iE2)' - (w2Si)2 



2] 



(100) 



Now consider the bulk terms (the second fine in (|47l) ) which are due to non-Hnear field interactions dilaton- two-form 
and dilaton-graviton. Their contribution is illustrated by the graphs / and 671, the second being pairwise. In quantum 
theory terms they can be interpreted as the coalescence of two virtual axions into the on-shell dilaton, and the 
coalescence of the off-shell dilaton and graviton into the on-shell dilaton. We have to compute the Fourrier-transform 
of the bulk current 



b 

J{x) 



247r 



H au\H' 



1 



da 



1 



(101) 



1 1 1 1 

Here in the products Hfu^xH'^" and (/yhfji, we have to substitute one field generated by the first string and another 
by the second one. 

In spite of the different different appearance of the bulk terms as compared to the world-sheet terms, one can cast 
them into the unique form |97|) as well. Consider the corresponding transformations for the first term in l|10ip depicted 
by the graph / on Fig. 4. The first order quantities to be substituted here read: 



2iTT 



'-^nq ''n ,4 



d'^q. 



(102) 



(We can omit the e-term in the denominator indicating the position of the pole since the resulting integral does not 
depend on its shift from the real axis). Now let us calculate the contraction 



-(27r/)2 



Ai(g)A2(g')g{^^^i.A}g'<'^V'2'^^^e'('=-'-'')%4 

d q 



q2qi2 



Ai(g)A2(fc - q)q{^Vi.x}{k - q){t^V^'^^}e'ik~iW)- 

q'^{k — g)2 ^ 



(103) 
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Phc. 5: The second order amplitudes of the two-form emission. The graphs a, b, c show contributions from deformations of the 
string world-sheets, the graphs d, e, f,h — the quadratic bulk terms. 



We are interested in the on-shell value of k'^, i.e. fc^ = 0. In this case one can write 

q^{k-qY \{q-kf J 2kq 



(104) 



Using this relation one can cast the above contraction into the form l|97|) with convention that the argument q'^ is 
used for the n—1 terms, and {k — q)^^ for the n — 2 terms. 

Performing similar calculations for the graphs (71,(72 and combining with the above, we obtain the total bulk term 

, , b 

in the form J93) with 6 ==9: 



(Yik) 
[kq) 



(uiU2)(SlS2) -I- (EiM2)(uiE2) 



2 ^{kU2f - {k^2? + {ku,f - {k^,Y 



(kq) 



(105) 



b b 

with the rule of getting 02 from Qi the same as before. In obtaining this expression we have used the delta-functions 
in the integrand fixing (qui) = (qTii) = and (9^2) = {ku2), (q'^i) = (fcSi). 



B. Two- form 



The source term in the equation ^ for the second order two-form field can be presented as the sum of ten graphs 
shown on Fig. 5. Here we have string contributions corresponding to exchange by the dilaton, two-form and graviton, 
the bulk terms describing coalescence of virtual quanta to the on-shell axion, but no contact term because of absence 
of multi-leg vertices associated with two-form. The string terms a„,6„,c„ are computed as the Fourrier-transform of 
the currents 



St''" 

Jn ix)^f 



Oil, 



(106) 



where perturbations of the mapping functions X!^, n = 1,2 are generated by the partner string n = 2,1 respectively. 
Using the same rearrangements as before, we obtain 



(107) 



where 



er(9) = 



[((?Ei)2 - (qu^y 



(108) 
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Phc. 6: The second order amplitudes of the graviton emission. The sum of all graphs is zero on the graviton mass-shell. 



and to get the second string term we must interchange indices 1^2 and momenta g'^ ^ fc^ — q'' . 
The bulk terms (i„,e„ are the Fourrier-transform of the bulk current 



(109) 



where again we have to take the products of fields generated by different strings. The result can be cast into the form 
(flOTl) with 



er 



SirffiG 



Vrikq) + {q^''ui''^){j:iq) - (gl^Si'^l ) (uiq) - Vr {kq) + (9^2"^] ) (Eafc) - (g[^S2'^l)(u2fc))' 



vrikq)) - (g[^Ui''l)(Sig) + (g[^Ei'^l)(^.ig) - V^ikq) - {q^^U2-^){^2k) + {q^n^2''\u2k)] , (110) 



2{kq) 

and the same rule for the second term. 



Graviton 



The source l|52|) in the second order graviton equation can be treated along the same line. It includes contributions 
of the thirteen graphs shown on Fig. 6. It turns out, however, that the projection of the total on-shell (fc^ = 0) 
current onto the graviton transverse polarization states gives zero. The reason lies in the dimensionality of the space 
transversal to the string: the configuration of the parallel strings is reduced to that of point particles in the 1+2 theory, 
where there is no on-shell transverse gravitational degrees of freedom. As it was shown in [s^l, the transformation 
to the parallel string configuration is always possible for a super luminal moving intersecting strings. So there is no 
gravitational radiation in this case and we do not give here the details of the calculation. 



VII. RADIATION 



Total radiation four-momentum loss can be presented in a standard way through the on-shell Fourrier transform of 
the source current, the vector k^ = (w, k), uj = |k| playing the role of the radiation four-momentum. In the case of 
the dilaton one obtains the following explicitly Lorentz-covariant expression: 



and similarly for the two-form 



IK/" h'^ 

% - ij k'^jj^^\Ji,m'sik^)d% 

P(B) - U k'^^^Mk)\'sik')d^k. 



(Ill) 



(112) 
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Alternatively, the latter quantity can bo presented as a square of the polarization projection of the current. Indeed, 
in three space dimensions the two-form field propagating along the wave vector k has a unique polarization state 



1 



6i 



ij = ^{e'(e^-efe'j), i,i = 1,2,3, (113) 

where and e*^ are two unit vectors orthogonal to k and to each other: 

e'^-e'' = 0, k-e^=k-e^ = 0. (114) 
Using antisymmetry and transversality of the two-form current fc'' J^,^(fc) = 0, and the completeness condition 

e^ef + ef ej = 6ij - kikj/oj\ (115) 

one finds 

Pi% = \j fc''|^|j;^)(fc)e,,f 5(fc2)d4fc. (116) 
Integrating over A;°, we finally obtain 

^r^) = -j piJ(^)(k)i^rf3fc, (117) 

^r^) = kj ^\\JiBm'd^k, (118) 

where 

J(B){^) = J'i){k)eii (119) 

with A:° = |k|. In what follows we shall use the parametrization of three- vectors by the spherical angles: k = 
a;[sin^cos</3, sin ^ sin cos ^] , e''' = [— sin t^, cosc^, 0],e^ = [cos^cos<p, cos^siny>, — sin^]. 

A. Cerenkov condition 

With our conventions, the radiation amplitudes associated with the first string contain the integral 

j 5{qu2)5{qi:^)5[{k - q)u,]5[{k - q)Y.,] fiq^'^^ A^^ (^20) 

where = dt^ — d^, and f{q) - is some regular function of q. In the Lorentz frame where = (1,0,0,0), S^* = 
(0,0, 0, 1), = 7 (1> Oj — 'f^cosa, usina), S2 = (0, 0, sin a, cos a), one can integrate over q^ using 

S[{k-q)ui]=6{q° -w). (121) 

Two other (5-functions 

6{qu2) = j^^S{uj + qyV cosa — q^v sina), (122) 

(5(gE2) = 7^^(5(qy sina -|- q^wcosa), (123) 



can be used to fix the values of qy,qz- 



wcosa (jsma 
qy = , qz = • 124 

V V 



The remaining (5-function 



S[{k - = 5ikz -qz) = S (^kz - (125) 
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no more depends on q, but rather restricts the value of the wave vector of radiation kzi 

uj sin a 

kz = ^ = Vp (126) 

This is Cerenkov's condition for an emitted massless wave. Indeed, in our frame an effective source of radiation is 
moving along z-axis with the velocity Vp, thus the quantity 

cos 61= — = — (127) 

UJ Vp 

defines Cerenkov's angle of emission if > 1, i.e. the source is super- luminal. Thus, radiation arises if the string 
relative velocity v and the inclination angle a satisfy the Cerenkov condition 

sina<'i;. (128) 

Given v, this condition will always be satisfied for sufficiently small a. In particular, it is identically fulfilled for a = 0, 
i.e. for parallel strings. Moreover, it can be shown that, if Vp > 1, the Lorentz frame and the world sheet coordinates 
T„,cr„ always exist such that two string look parallel [37||. In the frame where the first string is at rest and streched 
along the z-axis, an effective source moves in z-direction, and radiation will be emitted along the Cerenkov cone 
around the z-axis with an angle 

V 

^'ch = arccos . (129) 

sma 

The remaining integral over qx in (|120p can be evaluated using contour integration. With our choice of coordinates, 
the scalar product g^d^ = —qxd, so we have the integral over q^ 



qi+p^ p 



where 

p=^q^^ + q2-q2 = ^^. (131) 

For the second term in Ip7|) containing the pole {k — q)^^ one obtains the same fixed values of q^,qy,qz, but the 
integration over qx gives the value 

qx ~ kx ^ -, = cosa + wsin6'sin</). (132) 

iv 

Summarizing the above results, we obtain: 

/ (^J^ + ^(k -~qf ) = [Ei^M + {lO'^E2Q2{k - 92)] 5{cose^ cos^ch), (133) 

where 

El = (S^di-'^d/(v-t) ^ j^^ _ ^ikd2-ujd^/v ^ (134) 

and the following complex vectors are introduced: 

qi'^ — — [w, —1/7, — cos a, sin a], (135) 

V 

k^-qt^ - —[0, i, 1, 0]. (136) 

V 

The presence of the delta-function at the right hand side of (|133p means that the total radiation loss is infinite. This 
could be expected since we deal with the stationary motion of an infinitely long strings. So it is natural to consider 
the radiation loss per unit length of the string at rest. Redefining the currents as 

J(k) = /(k)J(cos - cos 6ich), (137) 
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and using the identity 



(5^(cos6'- cos6lch) = ^(5(cos6'- cos6'ch), (138) 

ZTT 



where L is the normalization length, we find 

T'f^) = ^"'^(0) = ^ / fcn%(k)|'<^(cos0-cos0ch)rf'fc, (139) 
V^j^^ = L-^P^B) = l^ yfc^|/(B)(k)|25(cos0-cos0Ch)rf'fc, (140) 

B. Relativistic peaking and spectrum enhancement 

According to l|133p . there is a frequency cut-ofi' due to exponential factors (|134p . Actually the amplitude is 
superposition of two terms which can be associated with contributions of two strings (this can not be taken literally: 
in the second order of the perturbation theory the superposition principle does not hold) . Two terms exhibit different 
frequency cut-off. The first sting term has a cut-off 

^<^, (141) 

which does not depend on the radiation angle, while the second one exhibits a (p— dependent cut-off: 

0. < ^ = ^ . (142) 

d(cosQ; -)- wsin^sint/?) 

This means that the angular distribution of radiation on the Cerenkov cone is anisotropic. This feature becomes 
especially pronounced in the ultrarelativistic case. 
In view of the identity 

cos^ a = sin^ 6* + (143) 

which holds on the radiation cone, in the ultrarelativistic limit 7 — > cx) the quantity ^ has a sharp minimum at 
ip = — 7r/2 corresponding to the direction of the moving string in the rest frame of the first string: 

(l + /3Vcosa2)], (144) 

27^ cos a 

where {3 — it /2 + Lp <^ 1. Due to the factor in the second term in (|119p the Cerenkov radiation is peaked around 
the direction = — 7r/2 within the narrow angular region 

/3<7"^ (145) 

Moreover, the frequency range associated with the second terms is substantially larger in the ultrarelativistic limit 
than that associated with the first term. Indeed, if k = 7 cos a ^ 1, one has: 



so the frequency range extends up to the frequency 



w < ^ (147) 



in the angular region l|145p . Therefore, radiation exhibits relativistic peaking in the forward direction in the same way 
as radiation of the ultrarelativistic particle. This could be expected, since relativistic peaking has purely kinematical 
nature. 
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C. Dilaton radiation 



Collecting the world-sheet contributions to radiation amplitudes we obtain after integration over q: 



I (fc)+ / (k) 



a I a 



El 



-^'f E 



sm f 



E2 
7^ 



'jv sm f 



2^7 



sin + IK cos E2 ^vsinO + iKC 



sm ( 



7W sm ( 



Gtt ( Ei[{k sin (y5 — i cos 4>) + 71; sin6']- 



E2 



jv sm f 



(148) 
(149) 



Here the upper line corresponds to the dilaton exchange, the second line — to the two-form and the last line to the 
graviton exchange. The bulk contribution after integration reads: 



/ (fc) ~ — [Ei (7k^ — 1 — ijv sin 9 cos ip) + iE2jv sin 0e"^] 



7^ 



\E\ (71; sin + K sin — j cos 0) + £'2*6*"*] 27K'i; sin Q. 



(150) 



Consider first the Cerenkov's threshold v — sin a, when the radiation cone shrinks to ^ch — 0. This corresponds 
to K = 1, ^ = 1/7. Since in this limit kd — 0, the exponents become equal E\ — E2 and the bulk term vanishes. 
The graviton contribution to l|148p vanishes too, while the dilaton and two-form contributions differs only by the 
coefficients. The total radiation amplitude at the threshold will be given by 



/thr(fc) = ^Ei (3qV' - SttV') 



(151) 



Integrating the expression l|139p for ^ = (the energy loss rate) over the angles in Sk = uj'^dujd cos 9 d(p, we obtain 
the infrared-divergent integral over frequencies. Introducing the inverse correlation distance A as an infrared cutoff 
parameter one finds: 



In the BPS limit jTll) one has 



2 2 n 2.2N2 f°° doj Ujd 

^thr = -;-(3a - Stt / ) / — exp(-2 — ) 
47r JA-i ^ ^7 



Pthr = / — exp(-2— ) 

TT J^-1 UJ VJ 



(152) 



(153) 



Integrating this over frequencies we obtain the total radiation rate in terms of the integral exponential function (see 
Appendix): 



Ei 1 



2d 
v"fA 



For small impact parameters d ^ vjA this expression can be approximated by the logarithm 



\2d^J 



(154) 



(155) 



where C is the Euler constant, e'~^ — 1.781072418. For large impact parameters, d ^ vjA, radiation exponentially 
falls off: 



Ei 1 



2d \ 



2d 



' vfA ) vjA ^ \ vfA ) 



2d \ 



(156) 



For 6* 7^ the radiation amplitudes are more complicated. Note, that in our reference frame the first string is at 
rest while the second one is moving. For this reason the radiation amplitudes are not symmetric with respect to two 
strings. As it was observed in the previous section, in the ultrarelativistic case the frequency range associated with 
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the factor E2 is much larger than that that with Ei, and in this limit the E2 terms are dominant. Collecting the 
dominant terms we obtain for the total amplitude of the dilaton emission in the ultrarelativistic case 

/ k = + 27C - ttV' 1 + r-^ - z^72«sin0e^^ 

^■^t, Lo^ L 8 \7i;sm6' / \ wsmt^ 

+/x2G7ri7Dsin6'(l + 2K70e*'^|. (157) 

Substituting this into l|139p and taking into account the angular peaking near = — 7r/2 we obtain the spectral-angular 
distribution of radiation per unit length in the vicinity of this direction: 

exp , (158) 



where 



2 2 

Vli= AGi?k'^ + Tifn Vt2=2Gi?K^+TifK+'^-^. (159) 



Dividing by oj we will get also the number of emitted dilatons 

dM 1 dV 



diodfi Lu dudp 



(160) 



Note, that the graviton, axion and dilaton exchange terms entering into the above expressions for the parameters fli,2 
exhibit different behavior in the invariant Lorentz factor k, the dominant as k ^ 00 being the gravitational term. 
The spectrum exhibits an infrared divergence, and in the forward direction P — it extends up to w (^max, where 

7K 

^^max = -^- (161) 

Integrating over frequencies with the infrared cut-off we obtain the angular distribution of radiation: 

dV 2 2 ^2^' f3')\ .f, dil+P^JK) \ 

—— = 32a K — „ . — El 1, . (162) 

dp (1 + k2/32)4 y J V ; 

Since the integral exponential function decays exponentially for large argument, the total radiation is peaked around 
/3 = within the angle 

13 < Vtk- (163) 

One can also obtain the spectral distribution of radiation extending the integration domain over f3 in l|158p to the full 
axis in view of the exponential decay of the integrand: 

rfP ^ 27m^j^2^^(^) ^ 2n,n2F2{z) + nlF,{z)], (164) 

duj dio 07 

where z = ^ and the functions Fi are expressed in terms of the probability integral l|226p : 

Fi{z) = (8z2 - 16;z + 30)-^ - (8z2 - 12z + 18- — )erfc(V^), 

VTTZ Z 



F2{z) = (-8z2-80 + 6)-^ + (8z2 + 12z-6+-)erfc(\/i), (165) 

•v/TTZ Z 



F^iz) = (8z2 + 32z + 6)-^ - (8z2 + 36z + 18- -)erfc(Vi). 

VTTZ Z 



For small frequencies u -^"fn/d these functions grows as 



Fi{z)^^, F2(z)~^, F3(z)^| (166) 
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while for large z they exponentially decay: 

, s 48e-^ hie-"" lOSe"^ 

Fi{z)^^=, F2{z)^--=, F^{z)^-=. (167) 
VTTZ'^ lyixz^ Z\/ttz° 

To regularize the infrared divergence one has to introduce the cut-off length A. After integration over frequencies uj 
from to infinity we obtain the total dilaton radiation rate in the ultrarelativistic limit: 

V = [{^Ihiy) + 2n,il2f2{y) + ^IIM)] , (168) 

where three new functions are introduced 

/i(y) = 5/(y) + crfc(V^)) (^-f Qy^ + 18y + ^) - - |y + 23^ , (169) 

/2(y) = /(y) - erfc(vy) dy' + Qy' -6y-^] + f ^y' + ^ ^"l ' (^^^^ 



/3(y) = /(y) + erfc(V^) (^^2/3 ^ i8y2 ^ ^ 1 j _ g^s + + 11 



where 



and the function /(y) is expressed through the generalized hypergeometric function l|230p 



For small y these functions grow as 



while for large y they tend to zero in view of the asymptotic relation (|233p . 
Similar expressions can be obtained for the total number of dilatons: 



where three more functions are introduced 



tending to zero at infinity, and 



(171) 



(172) 



/(2/) = 12^/f2^^2fi i;||-y) -31n(4ye^). (173) 



/i(y)~151ni /2,3(y)~31ni (174) 

y y 



M = + 20ir!2^2(y) + VlITM\ (175) 
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^i(y) - crfc(V?) ( 42/2 _ I2y - 39 + ) - ^ (4y2 _ - 30) - 6/(y), (176) 



^2(y) = erfc(Vy) V - 12y - 9 + - j + ^ (4y2 + iQy + 6) - 2/(y), (177) 
^3(y) = erfc(Vy) (^4y2 + 36y + 9 + ^ j - (4^^ + 34y - 6) - 6/(y), (178) 



^2,3 ~- (179) 

y y 



for small y. 

As we have noted, in the ultrarelativistic limit the dominant contribution to radiation amplitude comes from the 
graviton exchange term. Leaving only this contribution we find for the radiation rate: 

T't^ = \'KG''o?n^K'g{y\ (180) 
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where 
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g{y) = 25/(2/) + erfc(V^) ( ^2/' - 30y^ + 1142;+ — 



8 , 94 

772/ - ^2/ + 131 



Similary, for the dilaton number: 

/ ■7^\ p-y 

g^y) = erfc(Vy)) (^+4y2 - GOy - 183 + — j - (4y2 _ 62y - 150) - 38/(y). 

The corresponding numerical curves are shown on the Figs l7|8l For small y one has: 

P(« = 200^GW«'ln(^^), 



(181) 

(182) 
(183) 

(184) 
(185) 



where we introduced the dimensionless dilaton coupling constant. These quantities rapidly grow with increasing 
Lorentz factor of the collision 7 (recall that k = 7 cos a where a is the angle between the strings). Thus, Cerenkov 
radiation is greatly enhanced for ultrarelativistic velocities. 



D. Two-form radiation 



After integration over the momentum g, the string contribution to the two-form radiation amplitude (|107p takes 
the form: 



f (fc) = 4\/2^3j3 ^ sine ~ e-'^) \ + 

fi'jvsmo \ 7 ? / 



El {i sin ip ~ K cos f) 



E2 
7^ 



271; sin \ 
+ 4:\/2Tr^ f fiG"fv sin 9 i Ei{k cos ip + i sin ip) 



(kb "'f - ivj sin 9 - 2'y^£,v'^ sin^ 9 cos ip 

[i^v sin 9 + kb^"^ — 2 cos (^7^] 



E2 



The bulk amplitude l|107p will read: 



/ (fc) = ( ^^Il^yl +4^Tr^f ^Gjv sin 9 j (^^[{1 - 2-f^)-fv sin 9 cos ip - i] + ^{2v sin 9 cos ip - i^)' 



(186) 



(187) 
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Consider first the Cerenkov thresiiold v = sin a. In the above expressions the graviton exchange term in the string term 
vanishes, and the bulk term is zero. The remaining string amplitudes due to the dilaton and the two-form exchange 
simplify as follows: 

/,h,(fc) = -^^E, (SvrV' - 3a V') • (188) 

Substituting this into (|140p and integrating over the angles we obtain the spectral distribution of the two-form 
radiation on the Cerenkov threshold: 

2^2 







/ — exp 









^hr = ^ (SttV' - 3a V') / — exp ( -2^ ) , (189) 



where A the infrared cutoff frequency. After integration over frequencies one finds: 

^5 f6 / 2d 



For small impact parameters d <^ vjA this expression can be approximated by 



(190) 



for d 3> vjA one can use 



^ 4(2^)^6 /„^AX 
^thr«— :^ln(^^j, (191) 



^. / 2d \ 2d / 2d \ 

Now consider the case of the arbitrary Cerenkov angles in the ultrarelativistic limit. Leaving only the relativistic 
string contribution, we obtain: 

fE2 r4V2^V%,.2.,.,..„/, , 7r^a2 ^«;e-'^ + 7t>sin. 2 



-4%/27r^GM[i(l - K^) + K7i;sin6le"^) - ^7^^^]}. (193) 
Substituting this into (|140p we obtain the spectral-angular distribution in the vicinity of = 7r/2 

^^^3-3:3 = -. , .o^s>^4 exp — , (194) 



dLodl3 dujdf3 lu + \ 

where 

XD = -3-^, Xs = — , X/i = • (195) 



Integrating over the angles we get: 

dV Wir^d 
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[Fi{z)xl + F2{z){xhK^ - Xd)% (196) 



where 



dW 167r^d2f2 ,,,,,,,,,, , , 

[Fi{z)xl + F2{z)(xhn^ ~ xofl (197) 



F4 = (Sz^ + 8z + 6)-^ - (Bz^ + 122 + 6 - 3/z)erfc(y2), (198) 



and F2 is given by l|165p . For small and large frequencies one has: 
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respectively. Finally, integrating over frequencies we find the total two-form radiation rate and the number of axions: 

^(5) ^ IGttV! ^^y^^y^ ^ (^^^2 _ ^^)2;^(y)) ^ (200) 

= {^ITM + ixHn' XDfT,{y)) , (201) 



where 



For small y 



My) = f{y) + erfc(Vy) Qy^ + ey^ + 5y + - + T ^ + ^) ' ^^02) 

My) = erfc(Vy) (^4y2 + Uy - 3 + - {iy^ + lOy - 6) - 2/(y). (203) 



■^4--, /4(y)~31ni (204) 

y y 



while for large y both functions tend to zero. 

The result of the calculation of the axion Cerenkov radiation in the fiat space-time [i^l is reproduced putting 
XD ^ Xh = 0. It reads 



rr = Mv)- (205) 



In our case the dominant contribution comes form the graviton exchange. In the BPS limit one has 

8 , , 

XD = XB, Xh = rr^XB, (206) 

where a is the dimensionless dilaton coupling constant. For large k the leading term is 

' ^ 0-4 2 /2 0/ ■ 207 
In the most interesting case of small y f2 — fi, so the the ratio 



can be large, e.g. for a = 1 and k = 5 it is equal to 4 • 10''. 
The dominant term for the number of axions is 



(208) 



The numerical curves /2(y), •^^2(2/) are shown on Figs.[9l[l0l For small y one has: 



A^= o -4 2 -^^ly). (209) 



' = .-4. .2 111 -V I (210) 



= — rri^- 211 



VIII. COSMOLOGICAL ESTIMATES 



Cosmic strings are formed as a network of long strings of the size comparable to the horizon size. Colliding strings 
intercommute and form closed loop. At some stage the scale-invariant string network is formed consisting of long 
strings and loops which move freely with relativistic velocities. Evolution of cosmic superstring networks was recently 
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discussed in [20|-[2J|. In many respects cosmic superstrings are similar to the gauge theory cosmic strings, with some 
distinctions, however. In particular, for gauge theory strings the probability of the formation of loops P is of the 
order of unity, whereas for F-strings P ~ 10""^ and for and D-strings P ~ 10~^. The cosmic superstring network 
has a scaling solution and the characteristic scale is proportional to the square root of the reconnection probability. 
A typical separation between two long strings is comparable to the horizon size t (we use the standard cosmological 
units), C(i) ~ \/Pt. The results of numerical simulation show that the network of long strings reaches an energy 
density 



Ps 



P- 



Pt^ 



(212) 



Let us estimate the energy loss of long strings due to Cerenkov radiation of dilatons and axions. Consider an 
ensemble of randomly oriented straight strings moving chaotically in space. Let choose one target string between 
them and introduce the Lorentz frame where it is at rest. Other strings will have different orientations and velocities, 
and we can characterize them very roughly as moving in three orthogonal directions with equal probability. Since the 
dependence of the Cerenkov radiation on the incHnation angle a is smooth, we can use for an estimate the particular 
result obtained for parallel strings {a — 0) introducing an effective fraction ly of "almost" parallel strings taking into 
account the effect of the angular spread. Assuming N to be the number of strings in the normalization volume V — L^, 
we have to integrate the radiation energy released V in the collision with the impact parameter d — x over the plane 
perpendicular to the target string with the measure N/L^ ■ 2TTxdx. To estimate the radiation power per unit time we 
then have to divide the integrand by the impact parameter. Multiplying this quantity by the total number of strings 
N to get the radiation energy released per unit time within the normalization volume, we obtain for the Cerenkov 
luminosity: 



Vv—t;— 2TTdx. 

V 



(213) 



For BPS strings we use as V the leading relativistic terms (|180p and (|207p . Taking into account that the string number 
density is related to the energy density l|212p via 



V fiL' 

and assuming for a rough estimate L ^ A ^ t, where t is cosmological time, we obtain 



2/^3 -2 4 7 r» / \ ^ 



1 



where 



Siiw) 



9i.y)dy, S2{w) 



f2{y)dy, w = 
1 A 



(214) 

(215) 
(216) 

(217) 
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The exact values of these integrals are given in the Appendix and shown on Figs lll|12l Note that the realistic value of 
7 is of the order of unity, so applying our formulas obtained in the limit 7 1 is only an order of magnitude estimate. 

Now we can calculate the energy density ec of Cerenkov radiation as a function of time in the radiation dominated 
Universe. The energy density of massless fields scales with the Hubble parameter H as so we have to solve the 

equation 



dec 
IT 



= -AHec + Qc, 



where H 



Thus we obtain for the energy density of the Cerenkov dilaton and two-form radiation: 



'-C 



yTT G-'a /i 7 vSi{w) , 
-b2{w)- 



(218) 

(219) 
(220) 



where to - the initial time of the long string formation. Finally, using the approximate formulas i|236|237p valid in the 
relativistic case, we find: 



^ 8007r2G3aVVln7^^'^^ 



SB 



) 2i27r6/^7^1n7i^ ln(t/to) 



-4 ? 



Pt^ 



(221) 

(222) 



Comparing the last expression for axions with the fiat space result [40] we observe an enhancement due to higher 
power of the Lorentz factor and bigger numerical coefficient. This difference is due to the fact that here the main 
contribution to the first order interaction between the strings comes from gravitational force which is proportional to 
energy. In view of the previous analysis [i^ we can conclude that Cerenkov radiation from long strings is non-negligible 
effect in the cosmic superstring network. More detailed analysis will be given elsewhere. 



IX. CONCLUSION 



In this paper we have studied in detail Cerenkov radiation from moving straight strings interacting with dilaton, 
two-form and gravity. Formation of the faster-than-light sources in the system of randomly oriented moving straight 
strings is rather generic. As we have shown, these sources have a collective nature and arise due to deformations of the 
strings world-sheets caused by their interactions via massless fields. These deformations propagate with superluminal 
velocities if the inclination angle is sufficiently small, for parallel strings the source velocity is infinite. Radiation wave 
vectors lie on the Cerenkov cone in the same way as in the case of Cerenkov radiation of point charges in dielectric 
media. 
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One interesting feature related to dimensionality of a string compared to a point charge is the absence of gravitational 
radiation in four space-time dimensions. This is related to the fact that the space transverse to the straight string is two- 
dimensional, so the emitted massless fields must Hve in 1+2 dimensional space-time rather than in four-dimensional. 
As it is well-known, gravity in 1+2 dimensions does not contain free gravitons, this is why one can expect gravitational 
radiation from straight strings to vanish. In higher dimensions this objection does not work, so Cerenkov gravitational 
radiation can be expected in space-time dimensions higher than four. In four dimensions Cerenkov mechanism works 
for the dilaton and the two-form field which is equivalent to a pseudoscalar. 

To avoid complications due to possibility of "physical" string intersections (leading to the well studied processes 
of intercommutation and formation of loops) we consider the "collision" of strings moving in parallel planes. At 
each instant of time there exists a point of minimal separation between the strings, and it is this point which may 
propagate with the superluminal velocity. When interaction between the strings via dilaton, two-form and gravity is 
taken into account, strings get deformed in the vicinity of this point, these deformations contribute to an effective 
radiation source. Another contribution comes from tensions associated with the first order fields which give rise to 
these deformations. The string deformations give contributions locaHzed on the world-sheets, while the field stresses 
give bulk contributions. Both have the same order of magnitude. 

Cerenkov radiation from strings has some peculiar features in the highly relativistic region. We have shown that 
in this case radiation exhibits strong beaming on the Cherenkov cone in the direction of the fast string in the rest 
frame of the target string. The main radiation frequency is proportional to the inverse impact parameter, but in the 
ultrarelativistic case the spectrum is enhanced to high frequencies proportional to the square of the Lorentz-factor 
of the collision. It is shown that in this limit gravitational interaction between string dominates and gives the main 
contribution to the effective sources of dilaton and two-form radiation. 

Cerenkov's mechanism can be regarded as an analog of the bremsstrahlung of point charges in electrodynamics, 
which gives the main contribution to radiation in plasma. In the string case, however, there is another radiation 
mechanism due to existence of the internal string dynamics: radiation from oscillating loops. This effect is of the first 
order in couplings between the string and massless fields. Cerenkov radiation arises only in the second order in these 
coupHng, so presumably it is less important. But as we have shown here, it has stronger dependence on the Lorentz 
factor of the string collision, so it must become dominant for highly relativistic strings. Also, in the cosmic string 
network it is a pairwise effect which gives contribution to the radiation loss proportional to the square of the density 
of strings. Our rough cosmological estimates indicate that Cerenkov radiation is a non-neghgible effect in the cosmic 
string context indeed. 
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X. APPENDIX 



Here we collect properties of some special functions used in the main text. 



A. Integral exponential function 



The definition: 



oo 




(223) 



z 



The series expansion: 



Ei(l, z) = -C -\nz + z- z^l2 + z^ /Id, - 



(224) 



where C is Euler constant, e'^ = 1.781072418. 
The asymptotic expansion: 




(225) 



28 



The definition: 



The series expansion: 



The asymptotic expansion: 



B. Probability integral 



z 

2 f 2 

erf (z) = —= / e"* dt. 
Vi" J 



erf (z) = ^ (1 - zV3 + zVlO - . . .) 



erfc(z) = 1 — erf (z) 



In Sec. 7 we used the following indefinite integral: 



2z2 4z2 



113 3 



where 



2F2 (ai,a2;/3i,/32;a;) = ^ 



fc=0 



f3ik,f32k k\ 



2' 2' 



, Qifc = a(Q; + 1) • • • (a + fc), . 



(226) 



(227) 



(228) 



(229) 



(230) 



is the generalized hypergeometric function. To obtain an asymptotic behavior of the latter for z — > 00 we use the 
following identity: 



exp (-w) r/^\ ./^r^Al33 

Inw ^ dw = Inz crf(Vz) - 4W-2F2 -, -; -, -;-z 

WTTw y TT V 2 2 2 2 



which can be easily proved integrating by parts. Then taking into account the integral 

exp {—w) 



/ 



Ini 



-dw = -C- 21n2, 



/TTW 



we find that for z ^ 00 



(231) 



(232) 



(233) 



Integration of the functions g{y) defined in (|18ip and f2{y) defined in (|17Qp over y can be performed analytically: 

Siiw) = ^ giy)dy ^SOO^ 



, 1113 \ / 1 1 1 1 

2F2 -7:,-7:;77>7:;-"' -2^2 -77, -7:; 7:> 7:; 



2' 2'2'2' 



2' 2'2'2' 
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+ 75w (1 - In (iwe^)) + ( w + -w* + 57^^ - 10 

2 3 

1 e-y^ , ^ „ ,189 
^(8w3-124u;2 + 750w + 567) + — , 



189 



w 



)erfc(-^/w) 



12 ^ 

S2{w) = £ f2iy)dy = 12^ 



(234) 



,.4 , ?r,„..2 



F ( 1 111 

2-^2 I — — : 



2' 2'2'2' 



+ 3u; (1 - In (4we^)) + — (60w - 16u;'' + 72w^ ~ 48w^ - 9)erfc(V^) + 



1 e~^ Pu 3 
^ ^ (l6w^ + 40^2 - 84^ - 18) + - . 
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(235) 
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For small arguments the leading terms are 



Si{w) ~ -75w\-aw, (236) 
5*2 (w) ~ -Swlnw. (237) 
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